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PREFACE 


The  study  reported  herein  was  conducted  at  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  under  the  sponsorship  of  the  Office, 
Chief  of  Engineers,  Project  Uai6i102AT22 , "Theory  and  Principles  of  Re- 
inforced Earth,"  Task  A2,  Work  Unit  OOU . 

The  finite  element  and  parametric  studies  described  were  performed 
during  the  period  July  1975  through  February  19TT  by  Drs.  M.  M. 
Al-Hussaini  and  L.  D.  Johnson,  assisted  by  Mr.  Y.  S.  Jeng,  all  of  the 
Soil  Mechanics  Division  (SMD),  Soils  and  Pavements  Laboratory  (S&PL). 
This  report  was  written  by  Dr.  Al-Hussaini,  assisted  by  Dr.  Johnson, 
under  the  general  direction  of  Mr.  C.  L.  McAnear,  Chief,  SMD,  and 
Mr.  J.  P.  Sale,  Chief,  S&PL. 

Directors  of  WES  during  the  investigation  and  the  preparation  of. 
this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  J.  L.  Cannon,  CE.  Techni- 
cal Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Sl) 

UNITS  OF  MEASUREMENT 

U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (SI)  units  as  follows: 

Multiply To  Obtain 

inches 
feet 

square  inches 

square  feet 

pounds  (mass) 

pounds  (force)  per 
square  inch 

pounds  (force)  per 
square  foot 

kips  (force)  per 
square  inch 

pounds  (mass)  per 
cubic  foot 

inches  per  minute 

degrees  (angle) 


25.4  millimetres 

O.30U8  metres 

6.U516  square  centimetres 

0.0929030it  square  metres 

O.U535924  kilograms 

689^.757  pascals 

117.88026  pascals 

689U.757  kilopascals 

16.OI8U6  kilograms  per  cubic  metre 

25. h millimetres  per  minute 

O.OI7I+5329  radians 


FINITE  ELEMENT  ANALYEI„  OF  A REINFORCED  EARTH  WALL 


PART  I:  INTRODUCTION  # 

Background 

1.  Reinforced  earth  consists  of  a soil  mass  whose  engineering 
characteristics  and  performance  have  been  improved  by  the  introduction  of 
small  quantities  of  frictional  material  that  possess  a relatively  high 
tensile  strength.  Many  reinforced  earth  walls  are  designed  as  gravity 
structures,  the  stability  of  which  depends  on  the  frictional  stresses  be- 
tween the  reinforcing  elements  and  the  surrounding  soil  and  the  strength 
and  elastic  properties  of  the  reinforcing  material.  The  basic  design 
concept  of  a reinforced  earth  wall  assumes  that  the  soil  mass  is  in  ac- 
tive failure  governed  by  Rankine  earth  pressure  theory.  The  developed 
lateral  pressure  is  assumed  to  be  counterbalanced  by  the  friction  force 
that  develops  between  the  reinforcement  and  the  surrounding  soil. 

2.  An  instrumented  reinforced  earth  wall  was  constructed  at  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WE3)  to  detemine  the 
performance  of  such  a wall  during  construction  and  surcharge  loading.^ 

The  wall,  shown  in  Figure  1,  was  reinforced  with  galvanized  steel  strips, 
U in.*  wide,  0.02U  in.  thick,  and  10  ft  long,  spaced  at  intervals  of 

2 and  2.5  ft  in  the  vertical  and  horizontal  directions,  respectively. 
Three  galvsinized  strips  along  the  center  line  of  the  wall,  located  at 
elevations  1,  5,  and  9 ft  above  the  bottom  of  the  wall  (Figure  2),  were 
instrumented  with  full  SR-U  strain  gage  bridges  at  points  1,  2.5,  5^  and 
7.5  ft  from  the  face  of  the  wall.  At  the  elevation  of  each  instrumented 
tie  and  1 ft  away  from  the  face  of  the  wall  (i.e.,  skin  element),  two 
pressure  cells  were  placed  to  measure  the  induced  vertical  and  horizontal 
pressures  within  the  fill  material. 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Sl)  units  is  shown  on  page  U. 
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Figure  2.  Schematic  showing  the  location  of  the  electric 
strain  gage  bridges  along  the  center  line  of  the  retain- 
ing wall 


3.  The  skin  element,  which  comprised  the  exposed  surface  of  the 
wall,  was  made  of  Alcoa  Til  high-strength  aluminum  panels  normally  used 
in  the  expedient  construction  of  airfields.  Each  panel  was  2 ft  wide, 
12  ft  long,  and  1.6  in.  thick  and  can  be  connected  lengthwise  to  other 
panels  by  integral  hinge-type  connections.  Figures  1 and  2.  The  rein- 
forcing elements  were  fixed  to  the  skin  element  by  double  angle  con- 
nectors, Figure  1.  Further  details  of  these  components  are  discussed 
later.  A full  description  of  the  reinforced  earth  wall  is  presented  by 
Al-Hussaini  and  Perry. ^ 
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U.  The  reinforced  earth  structure  was  constructed  to  a height  of 
12  ft  and  loaded  by  a uniform  surcharge  pressure  to  failure.  Instrumen- 
tation readings  and  analyses  during  construction  and  surcharge  loading 
are  given  in  Reference  1.  The  exact  cause  and  pattern  of  failure  are 
not  known  since  failure  was  catastrophic;  however,  it  is  probable  that 
failure  was  initiated  either  by  fracture  of  one  of  the  connections  which 
joined  the  reinforcing  strips  to  the  skin  element,  by  failure  of  the 
skin  element  due  to  buckling  and  shear,  or  by  failure  of  the  reinforcing 
strips . 


Purpose  and  Scope 

5.  Accurate  analysis  of  a system  incorporating  reinforcing  strips 
placed  at  known  intervals  in  the  vertical  and  transverse  directions 
within  the  soil  mass  poses  a three-dimensional  (3D)  problem.  Proper 
modeling  of  a reinforced  earth  wall  also  involves  reasonable  simulation 
of  (a)  the  mechanical  behavior  of  the  earth  fill,  skin  element,  and  rein- 
forcing strips,  and  (b)  the  interaction  mechanisms  between  these  compo- 
nents of  the  wall.  Field  observations  during  and  u .bsequent  to  the 

construction  of  the  reinforced  earth  wall  at  WES  w re  not  in  complete 

2 

agreement  with  the  theoretical  method  proposed  originally  by  Vidal  and 
3 

later  by  Lee. 

6.  The  finite  element  (FE)  method  was  consequently  proposed  to 
simulate  as  closely  as  possible  various  sequences  of  construction  of  the 
reinforced  earth  wall  constructed  at  WSS  and  to  predict  the  distribution 
of  stresses  and  deformations  within  the  reinforced  earth  mass.  Although 
3D  formulations  and  codes  based  on  the  FE  method  are  possible,  their  use 
for  reinforced  earth  walls  would  involve  great  amounts  of  human  and  com- 
puter effort.  The  3D  problem  is  therefore  approximated  by  a structur- 
ally equivalent  two-dimensional  (2D)  system. 

7.  The  soil-structure  plane  strain  2D  FE  code  originally  devel- 

U 

oped  by  Clough  and  Duncan  for  analyses  of  Port  Allen  and  Old  River 
Locks  and  subsequently  modified  by  Radhakrishnan  and  Jones^  was  adopted 
for  use  in  this  study.  This  code  (see  Appendix  A),  which  uses  the 
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5-node  isoparametric  quadrilateral  element,  is  capable  of  simulation  of 
the  incremental  buildup  of  the  wall  and  the  concentrated  loading  con- 
struction sequences  needed  for  this  study.  The  nonlinear  behavior  of 
the  earth  fill  of  th ; reinforced  wall  is  simulated  by  a hyperbolic  for- 
mulation developed  by  Duncan  and  Chang. ^ The  behavior  of  the  metal  skin 
element  and  reinforcing  strips  is  assumed  linear  elastic  until  the  yield 
stress  is  reached;  thereafter,  these  metal  components  are  assumed  to 
fail  plastically. 

8.  Interface  elements  available  in  the  FE  code  were  also  used  to 
accommodate  slippage  between  and  separation  of  the  components  of  the  re- 
inforced wall.  The  interface  elements  are  especially  important  in 
separating  the  ends  of  the  reinforcing  strips  from  the  original  in  situ 
earth  face  which  supports  the  back  of  the  reinforced  earth  mass,  thereby 
permitting  the  stresses  in  the  reinforcing  strips  to  approach  zero  at 
the  ends. 

9.  FE  analyses  of  the  reinforced  wall  were  performed  for  the  con- 
struction of  the  wall  and  surcharge  loading  to  failure.  The  numerical 
predictions  of  the  FE  analyses  are  compared  with  field  observations,  and 
conclusions  regarding  design  analyses  are  presented. 
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PART  II:  MODELING  THE  MATERIAL  BEHAVIOR  OF 

THE  REINFORCED  EARTH  WALL 


10.  Because  of  the  composite  nature  of  the  reinforced  earth  wall, 
it  is  essential  that  the  stress-strain  relationships  for  the  foiondation 
soil,  ^ill  material,  reinforcing  strips,  ani  skin  element  are  accurately 
defined.  It  is  also  essential  that  the  interaction  between  these  com- 
posite materials  within  the  reinforced  mass  be  considered. 


Stress-Strain  Behavior  of  Foundation  Material 


11.  The  stress-strain  relationship  for  the  foundation  material 
was  approximated  by  the  hyperbolic  formulation,  summarized  in  Appendix  B, 
and  stated  by  the  following  expressions;* 

E^  = (1  - X^)^E.  (la) 


1 


1 


Rf(Ol  - - sin  (p) 

2c  cos  (}>  + 2o^  sin  <(> 


where 

E = tangent  modulus 

E.  = initial  modulus 
1 

K = hyperbolic  loading  parameter 
= atmospheric  pressure 
0^  = minor  principal  stress 
n = pure  number 
R^  = failure  ratio 


(lb) 


(Ic) 


* For  convenience,  symbols  and  abbreviations  used  herein  are  listed  and 
defined  in  the  Notation  (Appendix  C). 
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= major  principal  stress 
<{i  = angle  of  internal  friction 
c = cohesion 

12.  Laboratory  tests  were  not  conducted  on  the  foundation  material 
because  these  tests  were  not  required  in  the  conventional  design  analysis 
of  the  wall.  However,  several  triaxial  Q tests  previously  conducted  on 
soil  specimens  with  similar  characteristics  were  used  in  the  FE  analysis. 
The  stress-strain  curves  and  Mohr-Coulomb  criteria  of  the  ILL  soil  used 
in  the  FE  analysis  are  presented  in  Figures  3 and  , respectively.  The 
stress-strain  data  were  replotted  in  the  linear  form.  Figure  5,  which 
shows  that  most  of  the  data  points  of  each  test  fell  on  a straight  line, 
indicating  that  the  shape  of  the  stress-strain  curve  of  the  ML  material 
conforms  with  the  hyperbolic  representation.  From  Figure  5,  the  value 
of  the  initial  modulus  E^  , represented  by  the  inverse  of  the  intercept 
on  the  ordinate,  and  the  value  of  ( O-,  - o ) , represented  by  the 

inverse  of  the  slope  of  the  line,  were  obtained.  The  failure  ratio 
is  equal  to  j where  is  the  princi- 

pal stress  difference  at  failure,  as  illustrated  in  Figure  3.  A summary 
of  the  ML  soil  parameters  needed  for  the  hyperbolic  stress-strain  repre- 
sentation is  tabulated  below. 


"3 

psi 

E. 

1 

psi 

(°1  - l), 

psi 

psi 

«f 

(P  ) 

\ a/ 

6.95 

3,700 

27.78 

29.85 

0.93 

0.1*73 

20.81* 

5,550 

62.51 

70.1*2 

0.89 

1.1*18 

1*1.6? 

8,000 

111.12 

125.00 

0.89 

2.835 

69.1*5 

12,500 

168.75 

188.70 

0.89 

1*.790 

Average  0.90 


13.  The  values  of  the  loading  stiffness  K and  the  exponential 
constant  n needed  to  incorporate  the  effect  of  the  confining  pressure 
0^  in  the  hyperbolic  constitutive  model  are  obtained  by  plotting  E^ 
versus  (a_/P  ) as  shown  in  Figure  6.  The  best  fit  line  was  used  for  the 

J 8. 
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(cr,  -CT3),  PSI 


£|,  IN. /IN. 

Figure  3.  Stress-strain  curves  for  the  ML  soil 
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PSI 


10 
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Figvire  6.  Variation  of  the  initial  modulus  E.  versus  (a  /P  ) for  the 

ML  soil  ^ ^ ^ 


analysis.  Since  adequate  information  on  the  volumetric  behavior  of  the 
ML  soil  was  not  available,  Poisson's  ratio  v was  assumed  constant  at 
0.33.  Poisson's  ratio  at  failure  was  assumed  to  be  O-.hQ. 

Stress-Strain  Behavior  of  Sand  Fill 


► 

i 

1 

1 


iL.  The  fill  material  used  in  the  construction  of  the  reinforced 
earth  wall  was  clean  concrete  sand  with  particles  ranging  from  subangu- 
lar  to  angular.  More  detailed  physical  properties  of  the  sand  fill  used 
are  presented  in  an  earlier  report.^  The  same  procedure  of  hyperbolic 
stress-strain  analysis,  presented  previously  in  describing  the  behavior 
of  the  foundation  soil,  was  used  to  obtain  the  constitutive  equations 
for  the  sand  fill.  The  hyperbolic  paraneters  and  other  properties 
adopted  for  the  FE  analysis  were  as  follows: 


ii 


15 


r 


R. 


K 


n 

0.5 


’d 

pcf 


0.85  580  0.5  9T.2  0.30 

Poisson's  ratio  at  failure  was  assumed  to  be  0.48. 


degrees 

36 


c 

psf 

0 


cjj-.ress-Strain  Behavior  of  Reinforcinp;  Strips 


15.  The  steel  reinforcing  strips  used  in  the  earth  wall  were  made 
from  galvanized  zinc-coated,  24-gage  steel.  Each  strip  was  0.025  m. 
thick,  4 in.  wide,  and  10  ft  long.  The  stress-strain  curve  for  a 
1/2-in. -wide  strip  is  shown  in  Figure  J.  This  curve  shows  the  classical 


16 


1 


stress-strain  curve  with  modulus  of  elasticity  E = 31.1  x 10  psi  and 


yield  stress  o =51  ksi.  The  Poisson's  ratio 

y 

analysis  was  assumed  equal  to  0.28. 


needed  for  the  FE 


Elastic  Behavior  of  Skin  Element 


l6.  The  skin  element  was  made  of  Alcoa  Til  high-strength  aluminum 
panels,  2 ft  wide,  12  ft  long,  and  1.6  in.  thick,  connected  lengthwise 
to  each  other  by  a hinge  connector.  The  physical  and  elastic  properties 
of  the  panels  are  shown  below. 

2 

Moment  of  inertia  per  foot  of  panel  width  = 1 368  in. 

Section  modulus  per  foot  of  panel  width  = 1.396  in. 

Modulus  of  elasticity  ' ^ 


Weight  per  square  foot 
Poisson's  ratio 
Thickness  of  the  panel 
Thickness  of  the  sheet 
Yield  stress 


= 10  X 10  psi 
= U.6l  lb 
= 0.33 
= 1.6  in. 

= 0.1  in. 

= 35,000  psi 


IT.  The  reinforcing  elements  were  fixed  to  the  skin  element 
through  a connector  as  shown  in  Figure  8.  The  connector  consisted  of 
double  angles,  each  of  which  was  1-1/2  by  3 by  1/k  in.  and  12  in.  long, 
used  to  grip  the  reinforcing  strip  and  tie  it  to  the  skin  element  by 
two  l/lt-in.  bolts.  This  type  of  connection  was  simple  and  adequate  as 
long  as  the  skin  element  exhibited  very  small  or  no  bending  deformations. 

Behavior  of  Interface  Elements 


l8.  The  interface  element  formulation  introduced  by  Goodman  and 

7 

his  coworkers  was  incorporated  in  the  FE  program  used  in  this  study. 

The  constitutive  behavior  of  the  interface  elements  was  based  on  the 
hyperbolic  shear  stress-displacement  model  originally  formulated  by 
Duncan  and  Chang^  (Appendix  A)  and  summarized  by  the  following  equations: 
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2 C + a tan  6 
a n 


where 


K = tangent  shear  stiffness 
st 


K.  = initial  shear  stiffness 


K 

Y 


J 

w 


o 


n 

m 

T 


c 


dimensionless  number 
unit  weight  of  water 
normal  stress 
exponential  number 
shear  stress 
adhesion 


a 

6 = angle  of  friction 

The  normal  stiffness  of  the  interface  element  is  assumed  to  behave  elas- 

0 

tically  and,  given  a high  value  of  10  psf,  to  inhibit  mutual  penetra- 
tion of  the  materials  adjacent  to  the  interface  element.  The  tangent 
shear  and  normal  stiffnesses  were  each  set  to  a smai 1 residual  value  K 

r 

of  10  psf  if  the  normal  stress  of  the  interface  element  became  tensile, 

permitting  voids  to  develop  in  the  interface.  The  shear  stiffness  is 

also  reduced  to  the  residual  shear  stiffness  K if  the  shear  stress  at 

r 

the  interface  equals  or  exceeds  the  Mohr-Coulomb  strength. 


Interface  behavior 

between  reinforcing  strip  and  sand 

19-  The  nonlinear  shear  stress-displacement  relationship  between 
the  sand  and  the  reinforcing  was  determined,  prior  to  the  field  test, 
using  a specially  designed  shear  box.  The  testing  device,  shown  in  Fig- 
ure 9,  is  similar  to  the  direct  shear  box  except  that  the  lower  frame 
was  replaced  by  a sheet  of  the  galvanized  steel  glued  to  a wooden  block. 


SHEAR 

FORCE 


NORMAL  FORCE 


Figure  9.  Modified  shear  box  for  determining  skin  friction  angle 
between  the  sand  and  the  reinforcing  material 
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The  sand  was  compacted  to  a dry  density  of  101.7  pcf.  Three  direct 

shear  tests  were  conducted  with  normal  stresses  of  10,  50,  and  100  psi  , 

respectively,  at  a shear  deformation  rate  of  0.5  in.  per  min. 

20.  The  shear  stress-deformation  relationship  between  the  sand 

and  the  galvanized  steel  is  shown  in  Figure  10;  these  curves  provide  a 

reasonably  accurate  determination  of  the  initial  shear  stiffness  . 

The  ultimate  shear  stress  x was  obtained  from  the  inverse  of  the 

ult 

slope  of  the  linearized  shear  stress-deformation  relationship  as  shown 

in  Figure  11.  The  failure  ratio  was  obtained  by  dividing  the  peak 

shear  stress  t , as  measured  from  Figure  10,  by  the  ultimate  shear 
max 

stress  5 these  tests  the  value  was  approximately  equal  to 

unity.  The  angle  of  friction  6 and  the  adhesion  C were  obtained  by 

Q. 

plotting  the  shear  stress  at  failure  x versus  a as  shown  in  Fig- 

max  n 

ure  12.  A summary  of  the  parameters  needed  for  the  hyperbolic  shear 
stress-displacement  behavior  of  the  interface  elements  between  the  sand 
fill  and  the  galvanized  steel  is  presented  below. 
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21.  The  influence  of  stress  level  was  accounted  for  by  plotting 

the  initial  shear  stiffness  K.  versus  (a  /P  ) as  shown  in  Figure  13. 

1 n a 

The  experimental  data  fell  on  a straight  line  on  the  logaritlimic  scale, 
thus  enabling  accurate  determination  of  the  exponential  parameter  m 
and  the  dimensionless  number  defined  in  Equation  2b.  For  the  in- 

terface element  between  the  steel  strip  and  the  surrounding  sand,  m is 

equal  to  1.1  and  K is  equal  to  3280. 

J 

Interface  behavior 

between  skin  element  and  sand 

22.  No  direct  shear  test  was  conducted  to  examine  the  interface 
behavior  between  the  skin  element  and  the  sand  fill.  However,  the  non- 
linear parameter  obtained  previously  for  the  sand  on  galvanized  steel 
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SHEAR  DEFORMATION  A,  IN. 


Figure  10.  Shear  stress-deformation  relationship  between  the  sand  and 

the  galvanized  steel  strips 


.5 


5 


NORMAL  STRESS  EFFECT  Pg) 

Figure  13.  Variation  of  versus  (Oj^/Pg^)  for  sand  on 

galvanized  steel 


h 


was  used  to  describe  the  interface  behavior  between  the  skin  clement  and 


the  sand. 


Interface  behavior 

between  sand  and  natural  soil 


23.  The  nonlinear  constitutive  equation  presented  in  Equation  2 


was  used  to  describe  the  interface  behavior  between  the  sand  fill  and 


the  surrounding  natural  ground.  Three  direct  shear  tests  in  which  the 


lower  part  of  the  shear  box  was  filled  with  the  flL  soil  and  the  upper 


part  was  filled  with  compacted  sand  were  conducted,  and  the  linearized 
shear  stress-displacement  relationship  is  presented  in  Figure  ll+.  The 


figure  shows  that  the  experimental  points  deviate  slightly  from  the 
straight  line,  which  indicates  that  the  hyperbolic  representation  is 
only  an  approximation  for  the  behavior  of  this  material.  The  Mohr- 


Coulomb  failure  envelope  for  sand  on  silt,  presented  in  Figure  15, 
showed  an  angle  of  friction  of  3^  deg  and  a small  adhesion  of  0.05  psi . 
2U.  The  influence  of  the  stress  level  on  the  behavior  of  the  in- 


terface elements  between  the  sand  and  the  natural  ML  soil  was  accounted 


for  by  plotting  the  initial  shear  stiffness  K.  versus  (a  /P  ) as  shown 

1 n a 


in  Figure  l6.  The  figure  shows  that  the  experimental  points  deviate 
significantly  from  the  ideal  straight  line;  however,  the  best  fit  line 


was  used  for  evaluating  m and  K.  . A summary  of  the  parameters 

J 


needed  for  the  nonlinear  constitutive  equation  for  the  interface  element 


between  the  sand  and  the  natural  ML  soil  is  presented  below. 


degrees 


ll+,523 


25.  Comparison  between  the  shear  stress-deformation  curves  ob- 


tained by  experimental  tests  and  those  generated  by  the  hyperbolic  formu- 
lation, shown  in  Figure  17,  indicates  a reasonable  agreement  between  the 
two  respective  curves. 
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normal  stress  CTn.  PCI 


Figure  15.  Shear  stress  versus  normal  stress  at  failure  for 

sand  on  silt 
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Figure  l6.  Variation  of  initial  shear  stiffness  with  respect  to  (oj^/Pg) 
for  interface  element  between  the  sand  fill  and  the  natural  ML  soil 
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Figure  IT.  Shear  stress-deformation  relationship  along  plane  between 

the  sand  fill  and  the  ML  soil 


PART  III:  IDEALIZATION  OF  THE  REINFORCED  FARTH  WALL 


26.  The  reinforced  earth  wall  is  a complicated  3D  problem,  the 
complexity  of  which  is  compounded  by  the  introduction  of  the  steel 
strips  as  a reinforcing  element.  The  solution  of  such  a problem  re- 
quires a substantial  amount  of  human  effort  and  computer  time.  In  order 
to  arrive  at  a solution  with  reasonable  time  and  cost,  it  is  necessary 
to  introduce  a certain  amount  of  idealization  and  simplification  with 
respect  to  geometry  and  modeling  of  material  properties. 

Idealization  as  2D  Problem 


Modeling  of  reinforcing  strip 

27.  The  3D  reinforced  earth  wall  was  idealized  as  a 2D  plane 
strain  problem  by  assuming  that  the  reinforcing  strips  (see  Figure  1) 
are  replaced  by  a plate  extended  to  the  full  width  and  breadth  of  the 
wall.  These  plates  are  illustrated  in  mesh  I (Figure  I8)  and  mesh  II 
(Figure  19)  which  were  developed  for  this  study.  Assuming  that  the 
major  response  of  the  strips  is  provided  by  an  axial  stiffness,  the 
total  axial  stiffness  of  the  reinforcing  strips  S of  the  wall  is: 


A.E. 

j=i 


A E 
s s 


(3) 


where 

n = total  number  of  strips  in  each  row 

A = cross-sectional  area  of  the  reinforcing  strip 
s 

E^  = modulus  of  elasticity  of  the  galvanized  steel 
= length  of  the  reinforcing  strip 

28.  The  equivalent  stiffness  of  the  plate  that  substitutes  for 
each  row  of  steel  strips  may  be  defined  as 
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30 


Figure  l8.  FE  mesh  I employed  in  the  analysis  of  the  reinforced 

earth  wall 


where 


A = equivalent  cross-sectional  area  of  the  steel  strips 
e 

= equivalent  modulus  of  the  steel  strips 
= equivalent  length  of  the  steel  strips 
28.  Because  the  stiffness  S and  the  equivalent  stiffness 

should  be  the  same  and  also  because  the  length  and  the  equivalent 

length  remain  the  same,  thus  by  equating  Equations  3 and  , the 

equivalent  stiffness  E^  can  be  determined  as 


E 
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nA  E 
s s 


(5) 


29.  For  the  prototype  reinforced  earth  wall  the  cross-sectional 
area  of  the  reinforcing  strip  was  O.O96  sq  in.,  the  modulus  of  elastic- 
ity of  the  galvanized  steel  was  31.1  10^  psi,  and  the  number  of  strips 

in  each  row  was  6.  The  thickness  of  the  equivalent  reinforcing  strip 
elements,  as  shown  in  mesh  I (Figure  I8),  was  6 in.  and  as  shown  in 
mesh  II  (Figure  19),  0.02U  in.  The  width  of  the  equivalent  plate  for 
both  meshes  (i.e.  length  of  the  reinforced  earth  wall)  was  16  ft.  Sub- 
stituting the  appropriate  values  in  Equation  5,  the  magnitude  of  E^ 
can  be  determined. 


Mesh  I:  E = ~ ■ = 2.21+  x 10^  psf 

e 0.5  X 16 


Mesh  II:  E 


_ 6 X 0.096  X 31.1  X 10^  _ 


8 


0.002  X 16 


= 5.598  X 10  psf 


The  unit  weight  of  the  equivalent  reinforcing  strips  of  mesh  I was  as- 
sumed identical  with  that  of  the  sand  backfill  because  of  the  very  small 
cross-sectional  area  of  the  actual  reinforcing  strips.  The  unit  weight 
of  the  equivalent  reinforcing  strips  in  mesh  II  was  taken  as  1+92  pcf, 
which  is  the  unit  weight  of  the  actual  strips;  the  very  small  area  of 
these  strips  should  cause  little  additional  contribution  to  the  stresses 
from  gravity  using  mesh  II . 
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Modeling  of  skin  element 

30.  The  behavior  of  the  skin  element,  like  that  of  the  rest  of 
the  reinforced  earth  mass,  was  treated  as  a plane  strain  problem  whose 
major  response  is  provided  by  bending.  The  aluminum  panels  that  repre- 
sent the  skin  element,  as  shown  in  Figure  1,  were  replaced  by  an  equiva- 
lent beam  element  of  similar  deflection  response  depending  on  the  dimen- 
sions of  the  skin  element  in  the  FE  meshes.  Figures  l8  and  19.  The 
equation  used  to  satisfy  equal  bending  deformation  response  between  the 
prototype  and  the  model  was : 


where 

E = modulus  of  elasticity  of  the  aluminum  panel 

3. 

= moment  of  inertia  of  the  aluminum  panel  per  unit  width 

L = length  of  the  beam  (S  ) between  any  two  rows  of  reinforcing 
strips  (see  Figure  l) 

E = equivalent  modulus  of  elasticity 
e 

= equivalent  moment  of  inertia  per  unit  width 
= equivalent  length  of  the  beam 

31.  Since  the  beam  lengths  of  the  skin  element  and  its  model  are 
the  same  (i.e.,  L = L ),  the  equivalent  modulus  of  elasticity  for  the 

8L  0 

skin  element  can  be  expressed  by: 


E = E (7 

e a I 

e 

From  paragraph  l6,  the  modulus  of  elasticity  and  the  moment  of  inertia 
are  10  x 10^  psi  and  1.368  in.^,  respectively.  The  depth  of  the  equiva- 
lent beam,  as  shown  in  Figure  l8  for  mesh  I,  is  6 in.;  therefore,  the 
equivalent  moment  of  inertia  is  equal  to  216  in.  Substituting  the 

values  of  E , I , and  I in  Equation  7,  the  magnitude  of  E for 

3.  9.  0 0 

mesh  I can  be  obtained  as  follows : 


Eg  = = 6.33  X 10^  psi  (9.13  x lO^  psf) 


The  magnitude  of  E for  mesh  II  (Figure  19)  with  skin  element  thick- 

^6  8 
ness  of  1.6  in.  was  3.3^  x 10  psi  (U.8  x lo  psf). 

32.  The  aluminum  panels  of  the  skin  element  were  hinged  rather 
than  bonded  together  into  an  elastically  rigid  unit;  the  actual  modulus 
of  the  skin  element  may  consequently  have  been  substantiallpr  less  than 
that  assumed  for  the  FE  analyses  using  meshes  I and  II. 

33.  Two  iterations  were  made  for  computation  of  stresses  during 
each  construction  and  loading  sequence.  The  I’irst  iteration  indicated 
the  appropriate  modulus  (based  on  input  material  properties)  to  use  for 
computation  of  the  stresses  during  the  second  iteration. 

FE  Mesh  Design 

3^;.  In  addition  to  idealizing  the  reinforced  earth  wall  as  a 2D 

problem,  there  are  a number  of  other  factors  which  affect  the  design  of 

the  FE  mesh.  Some  of  these  factors  are: 

a.  Design  of  elements  required  to  obtain  a satisfactory 
result . 

* b.  Boundary  conditions. 

c_.  Wimiber  of  construction  and  loading  sequences. 

Location  of  interface  elements. 

Design  of  elements 

3(j.  It  is  always  desirable  to  use  a minimum  number  of  elements 
that  achieve  a satisfactory  degree  of  accuracy  so  that  the  computer  time 
and  hence  computation  cost  are  minimum.  However,  for  a given  number  of 
elements,  the  degree  of  accuracy  in  simulating  a particular  behavior  de- 
pends upon  the  correlation  between  the  displacement  pattern  assumed  and 
the  displacement  of  the  actual  events.  Because  of  the  uncertainties  re- 
garding the  assumed  displacement  and  the  preliminary  nature  of  this 
study,  simple  meshes  using  relativel:,  targe  numbers  of  rectangular  ele- 
ments, as  shown  in  Figures  l8  and  19,  were  developed  to  improve  accuracy. 

36.  The  elements  in  mesh  I,. Figure  I8,  were  square  to  minimize 
distortion  of  results  that  might  occur  if  the  aspect  ratio  is  different 
from  unity.  The  elements  in  mesh  II  were  reduced  in  total  number 
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compared  to  mesh  I to  increase  computation  efficiency  and  dimensions 
were  adjusted  to  closely  simulate  actual  dimensions  of  the  components  of 
the  wall. 

Boundary  conditions 

37.  The  boundaries  should  be  located  so  that  they  will  have  neg- 
ligible effect  on  the  main  structure  to  be  analyzed.  In  view  of  the  pre- 
liminary nature  of  this  analysis,  the  boundaries  were  kept  at  a distance 
of  2 ft  from  the  reinforced  earth  mass.  No  calculation  was  used  to  de- 
termine the  influence  of  the  boundary  conditions  on  the  FE  results. 

38.  The  left  and  right  vertical  boundaries  were  confined ■ later- 
ally and  the  bottom  horizontal  boundary  was  confined  both  laterally  and 
vertically.  During  each  construction  sequence  of  the  wall,  the  lateral 
displacement  was  set  to  zero  to  simulate  actual  construction  in  which 
the  lateral  displacements  of  the  skin  element  were  restricted. 
Construction  and  loading  sequences 

39-  Twelve  buildup  increments  1 ft  high  each  were  used  to  con- 
struct the  wall  described  by  mesh  I,  while  six  buildup  increments  2 ft 
high  were  used  to  construct  the  wall  described  by  mesh  II.  Concentrated 
loads  following  construction  were  placed  on  each  nodal  point  at  the  top 
©■f  the  wall  mesh  to  cause  surface  pressures  of  250  psf  per  increment. 

Six  increments  were  applied  to  reach  the  failure  condition  of  I5OO  psf 
observed  with  the  actual  wall. 

Location  of  interface  elements 

UO.  The  interface  elements  within  the  FE  meshes  are  presented  in 
Figures  I8  and  19-  These  interface  elements  increase  flexibility  in  ac- 
counting with  the  interaction  effects  between  the  different  components 
of  the  wall  such  as  slippage  and  separation  between  the  different  com- 
ponents. Mesh  I contains  vertical  interface  elements  between  the  skin 
element  and  backfill  material  and  between  the  back  face  of  the  wall  and 
the  original  in  situ  earth.  Horizontal  interface  elements  are  located 
between  the  reinforcing  strips  and  backfill  material  and  between  the 
backfill  and  original  ground  surface  of  the  wall.  Mesh  II  contains 
vertical  interface  elements  between  the  back  face  of  the  wall  and  the 


original  in  situ  earth  and  horizontal  interface  elements  between  the 
reinforcing  strips  and  backfill  material. 

Examination  of  the  FE  Mesh 


Ul.  Two  FE  analyses  were  performed  using  mesh  I and  mesh  II  in 
order  to  determine  the  influence  of  the  shape  of  the  FE  mesh  on  the 
results.  The  results  of  the  FE  analyses  obtained  from  both  meshes  were 
compared  and  correlated  with  field  measurements  at  the  end  of  construc- 
tion as  well  as  prior  to  failure.  End  of  construction  (EC)  refers  to 
completion  of  construction  of  the  reinforced  earth  wall;  prior  to  fail- 
ure (PF)  refers  to  placement  of  1500  psf  of  surcharge  load  on  the  top 
surface  of  the  wall. 


Tensile  stress  distribu- 
tion  in  reinforcing  strips 

h2.  The  tensile  stress  obtained  by  the  FE  method,  which  is  based 
on  a reinforcing  plate  covering  the  entire  area  at  a given  reinforced 
elevation  of  the  wall,  was  converted  to  tensile  stress  carried  by  the 
reinforcing  strips  as  actually  placed  in  the  field  by  the  following 
expression: 
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where 

a = tensile  stress  in  the  reinforcing  strip 

t^  = equivalent  thickness  of  the  reinforcing  tie  in  the  FE  mesh 
= equivalent  width  of  the  reinforcing  tie  in  the  FE  mesh 
N = number  of  reinforcing  ties  in  each  reinforced  elevation 
t = actual  thickness  of  the  reinforcing  tie 
W = actual  width  of  the  reinforcing  tie 

= equivalent  tensile  stress  computed  by  the  FE  analysis 
The  parameters  expressed  in  Equation  8 for  mesh  I and  mesh  II  are: 
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0.024  in 

4.00  in. 

6 

6.00  in. 

192.0  in. 


0.024  in. 
4.00  in. 

6 

0.024  in. 
192.0  in. 


43.  Comparisons  of  tensile  stresses  obtained  from  the  FE  method 
using  mesh  I and  mesh  II  and  those  obtained  from  actual  measurement  at 
the  end  of  construction  and  prior  to  failure  are  shown  in  Figures  20 
and  21,  respectively.  Figure  20  shows  little  difference  in  results  be- 
tween mesh  I and  mesh  II  at  the  end  of  construction.  The  tensile  stress 
distribution  from  the  FE  analysis  assumes  a parabolic  shape  similar  to 
that  obtained  from  actual  measurements  for  the  tie  at  elevation  1 ft 
above  the  base  of  the  wall.  The  locations  of  the  maximum  tensile  stress 
as  obtained  from  the  FE  analysis  agreed  well  with  those  measured  in  the 
field. 

44.  The  tensile  stress  distribution  along  the  reinforcing  ties 
prior  to  failure.  Figure  21,  indicates  that  the  FE  results  of  mesh  II 
are  slightly  higher  than  those  from  mesh  I.  However,  the  shapes  of  the 
tensile  stress  distribution  curves  for  mesh  II  are  much  closer  to  actual 
field  data  than  those  obtained  for  mesh  I.  These  results  indicate  that 
the  square  elements  as  represented  in  mesh  I reduce  the  maximum  tensile 
stress  as  well  as  the  curvature  of  the  tensile  stress  distribution 
curve. 

Lateral  pressure  distribution 

45.  The  lateral  pressure  distributions  along  a vertical  plane 
1.25  ft  from  the  skin  element  as  obtained  from  the  FE  analyses  of  mesh  I 
ard  mesh  II  are  shown  in  Figure  22.  These  results  indicate  that  while 
the  lateral  pressures  at  the  end  of  construction  are  slightly  higher  for 
mesh  I than  for  mesh  II,  the  trend  reverses  itself  prior  to  failure. 

The  resiolts  of  the  FE  analyses  for  both  EC  and  PF  cases  are  approxi- 
mately similar  but  somewhat  higher  than  those  obtained  from  the  Rankine 
active  earth  pressure  theory.  The  PF  lateral  pressure  near  the  top  of 
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Figure  21.  Influence  of  FE  mesh  on  tensile  stress  distribution  along 
the  instrumented  reinforcing  strips  prior  to  failure 
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b.  PRIOR  TO  FAILURE 


P’igure  22.  Influence  of  FE  mesh  on  lateral  pressure  acting  on 
vertical  plane  1.25  ft  from  skin  element 
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the  wall  as  obtained  from  the  FE  analysis  is  greatly  increased  for 
mesh  I,  a weakness  in  the  FE  program. 

U6.  From  the  foregoing,  it  appears  that  both  mesh  I and  mesh  II 
provided  approximately  the  same  results,  indicating  that  the  use  of 
square  elements  of  aspect  ratio  of  1.0  did  not  significantly  improve  the 
results.  Firrthermore , fewer  elements  are  contained  in  mesh  II,  sig- 
nificantly reducing  computation  and  increasing  accuracy.  Therefore, 
mesh  II  was  adopted  for  the  actual  FE  analysis  of  the  reinforced  earth 
wall . 
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PART  IV:  FE  ANALYSIS  OF  THE  REINFORCED  EARTH  WALL 


1*7.  Reliable  FE  analysis  of  tiie  reinforced  earth  wall  requires 
accurate  infomation  on  material  properties,  proper  modeling  of  deforma- 
tion and  failure  mechanisms,  and  proper  modeling  of  the  interaction 
effects  between  the  components  of  the  wall . The  previously  described 
models  of  material  behavior  (Part  II)  and  idealization  of  the  wall  by 
the  computer  code  (Part  III)  provide  only  a rough  approximation  of  the 
actual  field  situation.  The  results  of  this  analysis  and  comparisons 
with  field  observation  will  consequently  be  used  to  help  determine  more 
appropriate  parameters  for  material  properties  of  the  components  of  the 
wall  and  improved  models  that  may  better  represent  field  behavior. 

U8.  An  important  assumption  needed  to  obtain  the  FE  plane  strain 
solution  is  that  the  reinforcing  strips  cover  the  entire  length  of  the 
wall.  Such  an  assumption  increases  the  shear  resistance  between  the  re- 
inforcing element  and  the  surrounding  sand  by  a magnification  factor 
, where  m^  (see  Figure  l)  is  the  ratio  of  the  horizontal  cross- 
sectional  area  of  the  wall  to  the  total  area  of  the  reinforcing  strips 
at  that  elevation.  Thus, 
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Therefore,  the  increase  in  the  shear  stress  at  the  interface  element  of 
mesh  II  due  to  the  computer  code  is  compensated  for  by  reducing  the 
shear  stresses  in  Equation  2 by  m^  . Thus,  the  constitutive  equation 
for  the  interface  element  used  in  the  actual  FE  analysis  is: 
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Therefore,  the  material  parameters  input  in  the  code  are  = ^^10  and 
^ = 3 deg. 

Ii9.  It  must  be  noted  that  Equation  10  can  only  provide  a crude 
simplification  of  the  real  behavior  between  the  reinforcing  strips  and 
the  surrounding  soil.  More  study  is  needed  to  provide  a better  consti- 
tutive model  to  describe  the  behavior  of  the  interface  elements.  In 


K ^ is  set  to  a small  residual 

St 

value  K of  10  psf  whenever  the  shear  stress  at  the  interface  element 
r 


Equation  10  the  tangent  shear  stiffnes 

value  K of  10  psf  whenever  the  shea 
r 

exceeds  a Mohr-Coulomb  strength  envelope  with  cohesion  equal  to  C /m„ 

3.  I 

and  an  angle  of  friction  equal  to  6/m^  . 

50.  The  ribs  of  the  skin  element.  Figure  8,  are  extended  only 
along  the  length  of  the  wall;  thus,  they  contribute  very  little  to  the 
bending  resistance  of  the  skin  element  in  the  vertical  direction.  The 
moment  of  inertia  used  in  the  calculation  is  that  for  the  aluminum  sheet 
of  0.1-in.  thickness  rather  than  that  of  the  integral  aluminum  panel. 


FE  Analysis 

51.  Two  FE  analyses  were  performed  using  mesh  II.  In  the  first 

analysis,  case  , the  toe  of  the  wall  represented  by  point  A of  Fig- 
ure 2 was  assumed  fixed.  In  the  second  FE  analysis,  case  , point  A 

was  assumed  free  to  move  in  any  direction.  It  must  be  noted  that  under 
field  conditions  point  A is  neither  fixed  nor  free  and  the  true  answer 
may  be  between  the  two  extreme  cases.  The  lateral  pressures  exerted  by 
the  backfill  and  surcharge  loading  following  construction,  lateral  move- 
ment of  the  skin  element,  and  the  tensile  stress  and  deformation  be- 
havior of  the  reinforcing  strips  are  given  below.  Comparisons  are  made 
with  field  observations  for  the  EC  and  PF  cases. 

Distribution  of  lateral  pressure 

52.  The  distributions  of  lateral  pressure  along  a vertical  plane 

Ul* 


1 ft  from  the  wall  face  determined  from  the  FE  analysis  for  the  free- 
arid  fixed-end  cases  are  compared  with  field  observations.  Figure  23 
shows  the  variation  of  the  lateral  pressure  as  a function  of  the  height 
of  the  wall  at  the  end  of  construction.  The  differences  in  results  be- 
tween cases  and  at  the  end  of  construction  are  negligible  ex- 
cept at  the  base  of  the  wall  where  the  fixed  end,  case  , showed 

higher  lateral  pressure  than  the  free  end,  case  . The  results  of 

the  FE  analysis  at  the  end  of  construction  for  and  cases  are 

somewhat  higher  than  those  determined  by  Rankine  theory  and  those  ob- 
tained by  field  measurements. 


Figure  23.  Pressure  distribution  at  vertical  plane  1 ft  from 
skin  element  at  end  of  construction 


53.  Field  measurements  for  lateral  pressures.  Figure  23,  were 
obtained  directly  from  readings  of  earth  pressure  cells  placed  along  the 
center  line  of  the  wall  1 ft  from  the  skin  element,  and  also  from  those 
stresses  generated  in  the  instrumented  reinforcing  strips  within  a seg- 
ment of  the  wall  bounded  by  , Figure  1,  surrounding  the  strip. 

The  excessive  lateral  pressures  determined  by  the  FE  analysis  may  be  at- 
tributed to  deficiencies  in  the  computer  code  and  to  the  poor  modeling 
of  the  skin  element.  It  must  be  noted  that  the  panels  of  the  skin  ele- 
ment were  connected  together  by  horizontal  hinged  joints.  Figure  1, 
rather  than  continuous  bonded  joints  as  assumed  in  the  FE  analyses.  The 
crude  approximation  of  the  tangential  shear  stiffness  in  the  interface 
elements  between  the  reinforcing  strips  and  backfill  material  could  also 
have  contributed  to  the  difference  in  the  results. 

5^.  Figure  2^  shows  the  FE  results  for  lateral  pressures  along 
a vertical  plane  1 ft  from  the  skin  element  prior  to  failure.  The  re- 
sults indicate  that  the  lateral  pressures  obtained  for  the  fixed  end, 
case  , are  significantly  higher  than  those  obtained  for  the  free  end, 

case  Cg  . Both  analyses  deviated  significantly  from  the  linear  distri- 
bution of  Rankine  active  earth  pressure  theory  but  agreed  well  with  field 
measurement  at  the  upper  portion  of  the  reinforced  earth  wall.  It  is 
very  difficult  to  tell  which  FE  analysis  agrees  with  field  conditions  at 
the  middle  portion  of  the  wall  since  the  pressure  cell  reading  agreed 
very  well  with  case  while  lateral  pressure  calculated  from  the  ten- 
sile stress  of  the  instrumented  strips  agreed  with  case  . At  the 

lower  part  of  the  wall  it  is  clear  that  the  FE  analysis  of  case 
provides  better  agreement  with  field  measurement  than  does  that  of 
case  Cg  . 

55.  Significant  scatter  in  the  results  of  the  FE  analyses  was 
observed  from  one  vertical  plane  to  the  other;  the  scatter  became  ap- 
preciably worse  closer  to  the  skin  element.  The  scatter  in  results  is 
not  well  understood  but  was  probably  influenced  by  the  interaction  ef- 
fects between  the  reinforcing  strips,  skin  element,  and  the  backfill 
material. 
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Figure  2h.  Pressure  distribution  at  vertical  plane  1 ft 
from  skin  element  prior  to  failure 

Distribution  of  tensile 
strain  in  reinforcing  strips 

56.  The  reinforcing  strips  located  at  elevations  of  1 , 5>  and 
9 ft  along  the  center  line  of  the  wall  were  instrumented  with  full  SR-U 
bridges  on  both  surfaces  at  points  1,  2.5,  5,  and  7-5  ft  from  the  face 
of  the  skin  element.  The  measured  strains  at  the  top  and  bottom  of  each 
instrumented  point  were  averaged  to  eliminate  the  effect  of  bending 
strains;  the  average  strain  was  considered  equal  to  the  tensile  strain 
at  that  point  in  the  reinforcing  strip.  Only  the  strains  at  the  end  of 
construction  and  prior  to  failure  are  presented.  For  the  FE  analysis 


the  strains  were  calculated  from  the  average  deformations  of  the  nodal 
points  of  each  element  In  the  instrumented  strips.  It  has  been  noticed 
that  in  certain  cases  the  deformation  for  tlie  two  nodal  points  at  the 
lowei’  face  of  the  element  were  unbearably  excessive;  in  such  case,  only 
the  deformations  of  the  two  nodal  points  at  the  upper  face  of  the  ele- 
ment were  considered. 

57-  Comparisons  of  tensile  strains  from  the  FE  analyses  for  the 
free-end  and  fixed-end  cases  and  the  actual  strain  gage  measurements  of 
the  field  test  at  the  end  of  construction  are  presented  in  Figure  25. 

The  figure  indicates  that  the  FE  results  for  the  free-end  condition, 
case  Cg  , are  slightly  higher  than  those  for  the  fixed-end  condition, 
case  , but  much  closer  to  actual  field  measurements.  The  FE  analy- 

sis indicated  that  the  maximum  strain  in  each  reinforcing  strip  occurred 
at  points  closer  to  the  skin  element  than  those  obtained  from  actual 
measurement.  In  general,  the  maximum  tensile  strains  at  the  end  of  con- 
struction for  case  are  in  better  agreement  with  field  results  than 

those  obtained  from  case  . 

58.  The  trend  of  the  results  at  the  end  of  construction  differed 
significantly  from  those  obtained  prior  to  failure  of  the  wall.  Strains 
determined  from  the  FE  analysis  prior  to  failure,  shown  in  Figure  26, 
indicated  a maximum  strain  value  close  to  the  skin  element  that  de- 
creased gradually  until  it  vanished  at  the  free  end  of  the  reinforcing 
strips.  The  results,  unlike  the  end  of  construction  condition,  shoved 
higher  strains  for  case  than  for  case  . For  the  strips  at  ele- 
vation 5 ft  above  the  base,  both  and  cases  showed  similar 

values  and  their  maximum  values  were  in  good  agreement  with  that  ob- 
tained from  field  measurement. 

Distribution  of  tensile 
stress  in  reinforcing  strips 

59-  The  in  situ  tensile  stress  distribution  along  the  reinforcing 
strips  of  the  reinforced  earth  wall  was  calculated  from  the  strain  gage 
readings  and  the  modulus  of  elasticity  of  the  galvanized  steel  using 
modulus  of  elasticity  F,  = 31-1  10  psi.  Therefore,  the  variation  in 

the  field  tensile  stress  distribution  is  similar  to  the  variation  in 
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Figure  26.  Tensile  strain  distribution  along  instrumented 
ties  prior  to  failure 
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tensile  strain  distribution  as  long  as  the  stresses  in  the  reinforcing 
strips  are  within  the  elastic  range.  The  tensile  stress  at  any  point 
for  the  FE  analysis  was  obtained  from  the  equivalent  tensile  stresses 
using  Equation  8.  Comparisons  of  the  tensile  stress  distribution  along 
the  instrumented  strips  at  the  end  of  construction  between  the  FE  analy- 
sis and  field  test  results  are  presented  in  Figure  27-  The  figure  shows 
that  the  FE  results  for  the  free-end  case  and  the  fixed-end  case  are 
identical,  but  somewhat  different  from  field  measurements.  The  maximum 
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Figure  27-  Tensile  stress  distribution  along  instrumented 
ties  at  end  of  construction 


tensile  stress  for  the  FE  analysis  occurred  at  a point  much  closer  to 
the  skin  element  than  that  observed  in  the  field.  Both  the  Fp’.  results 
and  the  field  data  show  a gradual  decrease  in  the  tensile  stress  from  the 
maximum  until  it  reached  zero  at  the  free  end  of  the  reinforcing  strip. 

6o.  Comparisons  of  results  between  the  tensile  stress  prior  to 
failure  for  the  FE  analysis  and  actual  field  observations  are  presented 
in  Figure  28.  The  figure  shows  that  the  tensile  stresses  obtained  from 
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Figure  28.  Tensile  stress  distribution  along  instrumented 
ties  prior  to  failure 


thf  FF.  iiriitlysis  art-  (,’enerally  lower-  th-in  those  obtained  from  field  mea- 
surements and  the  rer-.ult.-  for  tire  fixed-erid  condition,  case  , are 

hiidier  than  those  for  the  free-end  condition,  case  • FE  analysis 

results  also  show  that  the  maximum  tensile  stresses  occurred  at  a point 
much  closer  to  tlie  ;-.kin  element  than  those  obtained  from  field  measure- 
ment. Reasons  for  the  lifferences  between  the  FE  results  and  field  ob- 
servations are  not  knowri  witfi  certainty;  however,  deficiencies  in  the 
computer  code  and  modeling,  difficulty  in  determining;  tire  actual  stiff- 
ness of  the  skin  element,  and  lack  of  precise  constitutive  equations  for 
describinf-;  behavior  for  the  materials  are  considered  major  contributing 
factors. 

Lateral  movement  of  skin  element 

61 . The  reinforced  earth  wall  was  constructed  such  that  each  sec- 

tion of  aluminum  paneling  added  to  fonri  the  skin  element  was  held  in 
place  to  prevent  lateral  movement  while  the  corresponding  layer  of  back- 
fill was  added.  Field  observations  indicated  negligible  deformations  at 
the  end  of  construction,  but  the  skin  element  deformations  became  sub- 
stantial prior  to  failure,  as  shown  in  Figure  29,  with  maximum  deforma- 
tion of  3.^  in.  at  the  top  of  the  wall.  Comparison  between  the  FE  analy- 
ses and  the  field  observations,  based  on  the  accumulative  deformation 
from  one  panel  to  the  next , shows  good  agreement  between  the  field  ob- 
servations and  FE  results  for  the  fixed-end  condition,  case  C.,  . For 
the  free-end  condition,  case  , the  FE  analysis  showed  excellent 

agreement  with  field  observations  at  the  lower  half  of  the  wall;  however, 
at  the  upper  half  of  the  wall  the  FE  resul-ts  deviated  significantly  from 
actual  observations  and  indicated  compression. 

Concluding  Remarks  on  the  FE  Analysis 

62.  Because  of  the  preliminary  nature  of  this  study  and  the  com- 
plex interaction  between  the  various  components  of  the  reinforced  earth 
wall,  it  i.s  considered  unlikely  that  total  agreement  between  the  FE 
analyses  and  field  measurements  could  be  obtained.  However,  it  is  shown 
that  the  boundary  conditions  of  the  wall  play  a significant  role  in  the 
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wall  deformation,  in. 


Figure  29-  Accumulated  lateral  skin  element  deformation 

prior  to  failure 


FE  results.  In  certain  aspects  such  as  the  deformation  of  the  skin 
element,  the  fixed-end  condition  showed  better  agreement  with  field  ob- 
servations, while  in  other  aspects  such  as  the  tensile  strain  in  the 
reinforcing  strips,  the  free-end  condition  showed  better  agreement. 
Other  problems,  such  as  improving  the  computer  code,  obtaining  better 
constitutive  equations  for  modelling  material  properties,  and  better 
representation  of  the  skin  element  stiffness  used  in  the  calculation, 
need  further  development  before  a total  a,-i-eement  can  be  achieved. 
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PART  V:  PARAMR;TRIC  STUDY 


General  Ob.jectives 

63.  In  the  previous  part  of  tiiis  repoi-t  it  was  shown  that  there 
are  some  differences  between  results  of  the  FE  analysis  and  ^hose  ob- 
tained from  field  observation.  The  differences  were  attributed  to 
many  factors,  some  of  which  were  discussed  previously.  Variables  which 
may  influence  the  performance  of  a reinforced  earth  wall  are  also  nu- 
merous. The  Influence  of  the  following  variables  on  the  tensile  stress 
disti-ibution  along  the  instrumented  reinforcing  strii')s  was  evaluated 
using  the  FE  analysis:  (a)  friction  between  the  reinforcing  strips  and 

the  fill  material;  (b)  end  conditions  of  the  skin  element;  (c)  stiffness 
of  the  skin  element;  and  (d)  length  of  the  reinforcing  strips. 

Influence  of  friction  between 
reinforcing  strip  and  fill  material 

6^4.  Equation  2 represents  the  constitutive  equation  which  is 
used  to  describe  the  stress-deformation  relationship  between  the  gal- 
vanized steel  and  the  concrete  sand  used.  Analysis  of  the  test  results 
presented  in  Figures  12  and  13  indicated  that  6 = I8  deg  and 

K.  = 3280  . These  parameters  were  reduced  by  a magnification  factor 
J 

m^  = 8 when  used  in  the  FE  code  to  obtain  tensile  stress  distribution 
for  cases  and  , shown  in  Figures  27  and  28.  In  order  to  de- 

termine the  influence  of  friction  between  the  reinforcing  strips  and 
the  fill  material,  FE  runs  were  performed  using  the  same  paruneters 
employed  in  cases  and  except  that  6 and  K.  were  not  re- 

duced. The  tensile  stress  distributions  for  the  fixed-end  condition, 
case  , and  the  free-end  condition,  case  , are  presented  in 

Figures  30  and  31,  respectively.  Comparison  between  Figures  30  and  31 
for  the  unreduced  interface  friction  and  Figures  27  and  28  for  the  re- 
duced interface  friction  indicates  that  the  increase  in  the  interface 
friction  increases  the  tensile  stress  in  the  reinforcing  strips;  the 
increase  in  stress  is  much  higher  for  the  free-end  condition,  case  Cj^  , 
than  for  the  fixed-end  condition,  case  . Based  on  this  analysis  it 
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Figure  30.  Tensile  stress  distribution  along  instrumented 
reinforcing  strips  at  end  of  construction  and  prior  to 
failure  for  case  (fixed-end  condition) 
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Figure  31.  Tensile  stress  distribution  along  instrumented 
reinforcing  strips  at  end  of  construction  and  prior  to 
failure  for  case  (free-end  condition) 
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can  be  concluded  that  the  increase  in  the  interface  friction  between 
the  reinforcing  strips  and  the  surrounding  fill  increases  the  maximum 
tensile  stresses. 

Influence  of  end 
conditions  of  skin  element 

65.  Two  boundary  conditions  of  the  skin  element  were  investigated 

to  determine  their  influence  on  the  tensile  stress  distribution  within 
the  reinforcing  strips.  In  the  first  analysis  the  lower  point  of  the 
skin  element  was  assumed  fixed,  cases  and  , and  in  the  second 

analysis  the  lower  point  was  assumed  free,  cases  and  . Com- 
parison between  cases  and  Cg  , shown  in  Figure  28,  indicates  that 

the  fixed-end  condition  yielded  higher  tensile  stresses  in  the  rein- 
forcing strips  than  that  for  the  free-end  condition.  However,  compari- 
son between  cases  and  Cj^  , where  the  interface  friction  between 

the  reinforcing  strips  and  the  fill  material  was  not  reduced,  shows 
that  the  free-end  condition  induced  higher  tensile  stresses  than  the 
fixed-end  condition.  Therefore,  it  is  not  possible  to  draw  definite 
conclusions  on  the  basis  of  boundary  condition  of  the  skin  element  alone 
since  the  interface  stress  between  the  reinforcing  strips  and  the  sur- 
rounding fill  plays  a major  role  in  dictating  the  magnitude  and  distri- 
bution of  the  tensile  stresses. 

Influence  of  skin  element  stiffness 

66.  In  the  previous  FE  runs  the  stiffness  of  the  skin  element  was 
based  on  the  moment  of  inertia  I of  the  metal  sheet  at  the  back  of 

X 

the  aluminum  panel,  on  the  assumption  that  the  T ribs  (see  Figure  2), 

which  were  placed  in  the  horizontal  direction,  did  not  influence  the 

stiffness  of  the  skin  element.  Knowing  that  the  thickness  of  the  metal 

sheet  is  0.1  in.,  the  moment  of  inertia  I per  foot  of  the  skin  element 
3 ^ 

is  0.0022  in.  To  evaluate  the  influence  of  skin  element  stiffness  on 
the  tensile  stress  along  the  reinforcing  strips,  another  FF  analysis, 

case  C-  , was  conducted.  Case  Cj-  was  similar  to  case  in  every 

^ ^ 3 

respect  except  that  the  moment  of  inertia  in  case  was  1.368  in., 

which  is  the  moment  of  inertia  per  foot  of  panel  as  shown  in  para- 

graph  16,  rather  than  0.0022  in.  , which  is  the  I^  for  the  aliiminum 
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sheet  only.  The  tensile  stress  distribution  along  the  instrumented 
strips  at  the  end  of  construction  and  prior  to  failure  for  case 
is  shown  in  Figure  32.  Comparison  between  results  for  case  and 

those  of  case  , shown  in  Figure  31,  indicates  that  the  magnitudes 
of  the  maximum  tensile  stress  are  not  significantly  influenced  by  the 
stiffness  of  the  skin  element.  However,  maximum  tensile  stresses  for 
low  skin  element  stiffness,  case  C|^  , occurred  at  points  much  closer 
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Figure  32.  Tensile  stress  distribution  along  instrumented 
reinforcing  strips  at  end  of  construction  and  prior  to 
failure  for  case 
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to  the  face  of  the  wall  than  those  for  high  skin  element  stiffness, 
case  . ITierefore,  it  can  be  concluded  that  the  stiffness  of  the 

skin  element  slightly  changes  the  shape  of  the  tensile  stress  distribu- 
tion along  the  reinforcing  strip  but  not  the  magnitude  of  the  maximum 
stress . 


Influence  of  length 
of  reinforcing  strips 

6T-  The  influence  of  the  length  of  the  reinforcing  strips  on  the 
magnitude  and  distribution  of  tensile  stresses  was  examined  by  increas- 
ing the  length  of  the  reinforcing  strips  and  the  depth  of  the  wall  from 
10  ft  for  case  to  20  ft  for  case  Cg  (Figure  33).  All  other 

parameters  between  cases  and  Cg  were  identical.  Comparison  be- 

tween Figure  33,  case  Cg  , and  Figure  30,  case  , indicates  that  the 
increase  in  the  length  of  the  tie  influenced  the  distribution  but  not 
the  magnitude  of  the  maximum  tensile  stresses  along  the  reinforcing  ties 
at  the  end  of  construction.  However,  for  the  prior  to  failure  case, 
maximum  tensile  stress  along  the  reinforcing  strips  showed  a substantial 
decrease  with  increasing  length  of  the  strips.  Therefore,  it  can  be 
concluded  that  the  increase  in  the  length  of  the  reinforcing  strips 
tends  to  decrease  the  magnitude  of  tensile  stress,  and  thus  less  cross- 
sectional  area  of  strips  may  be  required  for  design  purposes. 


Concluding  Remarks  on  the  Parametric  Stud\ 


68.  It  has  been  shown  that  reinforced  earth  is  a complex  problem 
whose  complexity  arises  from  different  materials  with  different  proper- 
ties that  constitute  the  reinforced  earth  mass.  The  influence  of  only 
a few  variables  was  considered  in  the  parametric  study.  For  the  cases 
considered,  the  magnitudes  and  locations  of  maximum  tensile  stresses  are 
summarized  and  compared  with  field  data  in  the  following  tabulation.  It 
can  be  seen  that  none  of  the  FE  analyses  actually  duplicated  field  con- 
ditions. While  the  magnitude  of  the  tensile  stress  in  case  C^  is 
closest  to  the  field  observation,  the  location  of  the  maximum  stresses 
in  case  C^  correlates  most  favorably  with  the  actual  field  da^a. 
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DISTANCE  FROM  SKIN  ELEMENT,  FT 


Figure  33.  Tensile  stress  distribution  along  instrumented 
reinforcing  strips  at  end  of  construction  and  prior  to 
failure  for  case  Cg 


6l 


TENSILE  STRESS.  KSl 


Tie  Klgvation  AV>ve  the  Foundation 


1.0  ft 

5.0  ft 

9.0  ft 

Maximum 

Max imum 

Maximum 

Friction 
Angle 
6 , deg 

Strip 
Length 
L , ft 

Tensile 

Stress 

usi 

Location* 

ft 

Tensile 

Stress 

psi 

Location* 

ft 

Tensile 
Gt  resE 
psi 

Ix>cation* 

ft 

Field 

test 



10 

39,000 

2.00 

30,000 

1.50 

29,000 

3.50 

Case 

3 

10 

28,700 

0.73 

29,lil0 

1.50 

17,220 

2.50 

Case 

3 

10 

25,liO 

0.75 

26,160 

1.50 

6,770 

0.25 

Case 

^3 

18 

10 

31,0l0 

0.75 

31,980 

2.50 

26,260 

2.50 

Case 

18 

10 

37,020 

1.50 

11,280 

2.50 

23,650 

1.50 

Case 

18 

10 

37,600 

2.  50 

39.760 

3.50 

23,670 

2.50 

Case 

18 

20 

22,260 

1.50 

21,510 

5.00 

17,650 

5.00 

• The  location  of  the  maximum  tensile  stress  was  measured  relative  to  the  skin  element. 


PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


69.  An  FE  study  was  performed  on  a reinforced  earth  wail  to  in- 
vestigate the  feasibility  of  using  numerical  methods  in  predicting 
stresses  and  deformations  of  the  various  elements  of  the  wall.  The  wall 
used  in  the  FE  analyses  was  16  ft  long,  10  ft  deep,  and  12  ft  high;  the 
earth  backfill  was  reinforced  with  galvanized  steel  strips  1 in.  wide, 
0.025  in.  thick,  and  10  ft  long,  spaced  at  2-  and  2-5-ft  intervals  in 
the  vertical  and  horizontal  directions,  respectively.  Based  on  the  re- 
sults of  the  FE  analyses  of  this  study  and  results  collected  from  labora- 
tory and  field  tests,  the  following  conclusions  are  drawn: 

The  one-dimensional  interface  elements  vertically  in- 
serted between  the  foundation  soil  and  the  back  of  the 
reinforced  earth  wall  are  necessary  to  permit  separation 
between  the  reinforcing  strips  and  foundation  soil  such 
that  the  tensile  stresses  developed  in  the  reinforcing 
strips  decrease  to  zero  at  the  back  of  the  wall  as  ob- 
served from  field  measurements.  Interface  elements 
horizontally  inserted  between  the  reinforcing  strips  and 
backfill  of  the  wail  are  necessary  to  model  the  slippage 
and  friction  forces  between  the  reinforcing  strips  and 
backfill . 

b.  The  assumption  of  using  the  FE  plane  strain  analysis  to 
solve  the  3D  problem  of  a reinforced  earth  wall  by  in- 
creasing the  area  of  the  reinforcing  strip,  to  cover  the 
entire  length  of  the  wall,  and  decreasing  the  shear  re- 
sistance between  the  fill  and  the  reinforcement  by  the 
same  proportion  is  a practical  and  satisfactory 
assumption. 

£.  The  FE  analyses  of  both  fixed-end  and  free-end  skin  ele- 
ments are  necessary  to  obtain  good  correlation  with  field 
data.  It  is  observed  that  the  FE  results  of  tensile 
stresses  along  the  reinforcing  strips  and  the  lateral  de- 
formation of  the  skin  element  showed  better  agreement 
with  field  observation  for  the  fixed-end  condition  than 
for  the  free-end  condition.  However,  field  measurements 
of  strains  along  the  reinforcing  strips  agreed  better 
with  the  FE  results  of  the  free-end  condition  than  with 
those  obtained  from  the  fixed-end  condition.  While 
lateral  pressure  cell  measurements  agreed  better  with  FE 
results  of  the  fixed-end  rather  than  free-end  conditions. 
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the  opposite  was  observed  when  comparison  was  made  with 
lateral  earth  pressure  calculated  from  tensile  stresses 
in  the  reinforcing  strips. 

ji.  The  magnitude  and  distribution  of  tensile  stresses  along 
the  reinforcing  strips  increase  with  increasing  interface 
friction  between  the  strip  and  surrounding  fill  material, 
decrease  with  increasing  reinforcing  strip  length,  and 
show  no  significant  change  with  change  in  skin  element 
stiffness . 


Recommendations 


70.  Based  on  the  findings  of  this  study,  it  is  recommended  that 

the  study  be  expanded  as  follows : 

It  is  desirable  to  improve  the  computer  code  and  make  it 
more  accurate  and  efficient  by  introducing  elements  that 
better  simulate  actual  events.  The  variable  node  of  bend- 
ing element  is  recommended  because  fewer  elements  may  be 
needed  and  oscillation  observed  in  results  using  the 
isoparametric  elements  may  be  reduced  or  eliminated. 

b.  Better  constitutive  equations  that  not  only  simulate  the 
behavior  of  each  material  in  the  reinforced  earth  mass 
tut  also  describe  the  interaction  between  different  mate- 
rials should  be  developed. 

c_.  More  studies  are  needed  to  model  and  improve  the  behavior 
of  interface  elements  between  the  reinforcing  strips  and 
surrounding  soil. 

d.  Hinges  between  panels  that  form  the  skin  element  need  to 
be  considered  in  the  PE  code. 
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Homop;eneous  Soil  Element 
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1.  A simple  and  practical  expression  for  describing  the  nonlinear 

8* 

soil  behavior  has  been  suggested  by  Kondner  who  assumed  that  the 
stress-strain  curve  as  obtained  from  triaxial  compression  tests  is  the 
following  hyperbola: 


^<^1  - "3^  = rrn: 


(bi) 


where  (a^  - o^)  is  the  principal  stress  difference,  is  the  axial 

strain,  and  a and  b are  parameters  that  can  be  determined  by 
experiment . 

2.  The  physical  meaning  attached  to  a and  b can  be  seen  in 
Figure  Bl  in  which  a is  equal  to  the  reciprocal  of  the  initial  tangent 
modulus  and  b is  equal  to  the  asymptotic  value  of  the  stress 


Figure  Bl.  Kondner  hyperbola  for  representation  of 
stress-strain 


* Raised  numerals  refer  to  similarly  numbered  items  in  the  References 
at  the  end  of  the  main  text. 
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difference  called  the  ultimate  stress  difference  - 03^^^  . Equa- 
tion B1  can  be  simplified  by  expressing  c^/(o^  ~ as  a linear  func- 

tion of  , which  would  enable  direct  evaluation  of  the  parameters  a 
and  b as  depicted  in  Figure  B2.  The  linear  form  of  Equation  B1  may  be 
written  as: 


a + be^ 


(B2) 


3.  Duncan  and  Chang^  expanded  Kondner's  hyperbolic  stress-strain 
relationship  and  utilized  it  for  incremental  FE  analysis  by  expressing 
the  tangent  modulus  E^  as  a function  of  (a^  - o^)  , 

Mohr-Coulomb  strength  parameters  c and  4>  as: 


E 


t 


1 


R^(l  - sin  4))(a^  - 
2c  cos  (j)  + 2a^  sin  4) 


(B3) 


in  which  is  the  ratio  of  (a^  - a^)  at  failure  to  the  asymptotic 

value  of  the  stress  difference  ( '^l  - '^sjult  • 


Figure  B3.  Linear  representation  of  the 
hyperbolic  stress-strain  function 


Another  characteristic  of  many  soils  is  the  dependency  of  the 

9 

initial  tangent  modulus  E.  upon  the  confining  pressure  o,  . Janbu 

i 6 ^ 

has  suggested  and  Duncan  and  Chang  have  employed  the  following  exponen- 
tial relation  between  the  confining  pressure  and  : 

(BU) 


where  is  the  atmospheric  pressure  expressed  in  the  same  units  as 

, K is  a dimensionless  number,  and  n is  a dimensionless  exponent. 

5.  The  hyperbolic  stress-strain  relationship  as  expressed  by  Equa- 
tions B3  and  bU  is  useful  and  simple.  The  parameters  K , n , c , 4>  , 

and  are  all  readily  determined  by  conventional  triaxial  compression 

tests . 

Interface  Soil  Elements 


6.  Although  the  stress-displacement  curve  from  direct  shear  tests 

cannot  be  used  directly  for  determining  the  stress-strain  relationship 

for  soil,  it  seems  likely  that  a curve  for  such  a test  is  indicative  of 

the  form  of  the  stress-strain  relationship.  A comprehensive  study  by 

k 

Clough  and  Duncan  has  shown  a resemblance  of  the  shear  stress  relative 

g 

to  the  displacement  to  the  hyperbolic  form  suggested  by  Kondner,  which 
is  appropriate  for  characterization  of  the  behavior  of  interface 
elements. 

T.  Typical  shear  stress  t versus  displacement  for  a direct 

shear  test  is  presented  in  Figure  B3a,  where  K.  and  K , are  the 

1 St 

initial  and  tangent  shear  stiffnesses,  respectively.  The  hyperbolic 
form  for  the  shear  stress  as  a function  of  the  relative  displacement  is; 


T 


A 


s 


a + bA 

s 


(B5) 


where  a and  b are  parameters  that  can  be  obtained  from  laboratory 
tests.  Equation  B5  may  be  transformed  to  a linear  form  as  shown  in 
Figure  B3b: 


bU 


m 


— = a + bA 
T s 


(B6) 


8.  It  has  also  been  suggested  by  Clough  and  Duncan  that  it  is 
desirable  to  determine  the  straight  line  for  the  test  data  by  connecting 
points  on  the  curve  corresponding  to  TO  and  95  percent  of  the  shear  re- 
sistance mobilized.  The  constants  a and  b are  related  to  the  ini- 
tial shear  stiffness  K.  and  the  ultimate  shear  resistance  t as: 

1 ult 


b = 


K. 

1 


ult 


(B7) 

(b8) 


The  relationships  expressed  by  Equations  BT  and  b8  can  be  used  to  deter- 
mine and  "^^^t  transformed  linear  hyperbolic  plot  shown 

in  Figure  B3b  by  simply  evaluating  a and  b . 

9-  The  tangent  shear  stiffness  K ^ at  : 

St 

bolic  curve  may  be  obtained  by  differentiating  Equation  B5  as: 


9-  The  tangent  shear  stiffness  at  any  point  on  the  hyper- 


K 


dT 

dA 


st  dA  , . , . ,2 

(a  + bA  ) 
s 


(B9) 


From  Equation  B5  it  is  possible  to  obtain  an  expression  for  A^  by  re- 
arranging terms  to  read; 


ar 


1 - bx 


(BIO) 


By  substituting  a , b , and  A as  expressed  in  Equations  BT,  B8,  and 

s 

BIO,  respectively,  in  Equation  B9,  an  expression  for  can  be 

obtained: 


K = K./l  - — — 

"V  Vt/ 


(Bll) 


Experiments  have  shown  that  ''’yq-j.  overestimates  the  actual  value  of  the 
shear  stress  at  failure  x^  , and  the  ratio  between  two  shear  stresses 


b6 


does  not  vary  significantly  for  a given  soil.  By  introducing  the 
failure  ratio  into  Equation  Bll , the  following  may  be  obtained: 


K 


st 


(B12) 


10.  From  Mohr-Coulomb  failure  criteria,  presented  in  Figure  b4 , 


Figure  bU.  Variation  of  friction  resistance 
with  normal  stress 

the  shear  stress  at  failure  may  be  related  to  the  normal  stress  by 

the  following  expression: 


= C + a tan  6 (B13) 

fan 

where  C is  the  adhesion  and  6 is  the  angle  of  skin  friction.  There- 
a 

fore,  by  substituting  Equation  B13  in  Equation  B12,  the  following  is 
obtained: 


K 


st 


R^t  > 

C + a tan  6 
a n ; 


(BlU) 


11.  Laboratory  studies  by  Clough  and  Duncan  have  shown  that  the 

initial  shear  stiffness  K.  may  be  related  to  the  normal  stress  o 

1 n 


by  the  following  eApunential  relationship: 


K.  = K,y  I r, 

1 J wlP, 


(B15) 


where  y is  the  unit  weight  of  water,  P is  the  atmospheric  pressure 

W O’ 

and  K and  m are  parameters  that  can  be  determined  from  laboratory 
tests.  Equations  Blit  and  B15  define  a simplified  and  practical  relation 
ship  which  describes  the  nonlinear  shear  stress-displacement  relation- 
ship for  the  interface  elements. 
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Experimental  parameters 

Equivalent  cross-sectional  area  of  reinforcing  strips 

Cross-sectional  area  of  reinforcing  strip 

Cohesion 

Adhesion 

Modulus  of  elasticity  of  aluminum  panel 
Equivalent  modulus  of  elasticity 
Initial  modulus 

Modulus  of  elasticity  of  galvanized  steel  reinforcing 
strips 

Tangent  modulus 

Moment  of  inertia  of  aluminum  panel  per  unit  width 

Equivalent  moment  of  inertia  per  unit  width 

Moment  of  inertia  of  metal  sheet  at  back  of  aluminum 
panel  per  unit  width 

Hyperbolic  loading  parameter 

Initial  shear  stiffness 

Dimensionless  number 

Residual  shear  stiffness 

Tangent  shear  stiffness 

Length  of  beam  (S2)  between  any  two  rows  of  reinforcing 
strips 

Equivalent  length 

Length  of  reinforcing  strip 

Magnification  factor 

Pure  number;  also  total  number  of  strips  in  each  row 
Number  of  reinforcing  ties  in  each  reinforced  elevation 
Atmospheric  pressure 
Failure  ratio 

Total  axial  stiffness  of  reinforcing  strip 

Equivalent  stiffness  of  plate  substituting  for  reinforc- 
ing strips 

Actual  thickness  of  reinforcing  tie 

Equivalent  thickness  of  reinforcing  tie  in  FE  mesh 

Water  content  of  soil 
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w 
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te 


Oi  - 03 


(»1  - »3) 

’1  - 03) 


ult 

T 


max 


ult 


Actual  width  of  reinforcinf’;  tie 

Equivalent  width  of  reinforcing’;  tie  in  EE  mech 

Dry  density  of  soil 

Unit  weight  of  water 

Angle  of  skin  friction 

Shear  deformation  or  displacement 

Displacement 

Axial  strain 

Stress  levels;  see  Equations  Ic  and  2c,  respectively 
Poisson's  ratio 
Normal  stress 

Tensile  stress  in  reinforcing  strip 
Yield  stress 

Equivalent  tensile  stress  in  reinforcing  strip  computed 
by  FE  analysis 

Major  principal  stress 

Minor  principal  stress 

Principal  stress  difference 

Principal  stress  difference  at  failure 

Ultimate  principal  stress  difference 

Shear  stress 
Peak  shear  stress 
Ultimate  shear  stress 
Angle  of  internal  friction 
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